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ABSTRACT 

HD 259440 is a BOpe star that was proposed as the optical counterpart to the 7-ray source HESS 
J0632+057. Here we present optical spectra of HD 259440 acquired to investigate the stellar pa- 
rameters, the properties of the Be star disk, and evidence of binarity in this system. Emission from 
the Ha line shows evidence of a spiral density wave in the nearly edge-on disk. We find a best fit 
stellar effective temperature of 27500—30000 K and a log surface gravity of 3.75—4.0, although our 
fits are somewhat ambiguous due to scattered light from the circumstellar disk. We derive a mass 
of 13.2—19.0 Mq and a radius of 6.0—9.6 Rq. By fitting the spectral energy distribution, we find 
a distance between 1.1—1.7 kpc. We do not detect any significant radial velocity shifts in our data, 
ruling out orbital periods shorter than one month. If HD 259440 is a binary, it is likely a long period 
(> 100 d) system. 

Subject headings: binaries: general - stars: emission-line. Be - stars: fundamental parameters - stars: 
individual (HD 259440, HESS J0632-K057) 



1. INTRODUCTION 



HD 259440 (=MCW 443, BD -^05 1291), which hes 
between the edge of the Monoceros Loop supernova rem- 
nant and the star forming regi on of the Rose t te Ne bula, 
is classified as a B Ope star by Morgan et al. j 1955[) and 



a BNO nne star by Turner 



(197B). Gutierrez-Soto et al 



( 2007 1 assign it a projected rotational velocity, V sini, ot 
430 km s~^. Otherwise, no detailed optical study of this 
star has been performed. HD 259440 has been proposed 
as the optical count erpart to a possible new 7-ray binary, 
HESS J0632-I-057 ( jAliaronian et al.||2007[ ). 7-ray bina- 
ries are a subset of high-mass X-ray binaries which also 
emit radiation in the MeV-TeV range. The very high 
energy emission in such objects is expected to originate 
from an interaction between the massive star's wind and 
the compact object, although the physical process es that 
produ ce this radiation are under debate (see e.g. Dubus 
2006| . 

Observations taken with the HESS telescope array be- 
tween 2004 March and 2006 March revealed the TeV 



point source HESS J 0632-H057 ( [Aharonian et al.||2007[ ). 
Hinton et al. (2009) followed up this detection with a 
26 ks observation of the field using XMM-Newton on 
2007 September 17. XMMU J063259.3-I-054801, the 
brightest source in the field, was positionally coincident 
with the HESS observation and centered on HD 259440. 
The X-ray source showed a gradual decline in flux over 
the course of the observation. X-ray observations ob- 
tained with Swift showed flux variations o n a variety 
of timescales, from several days to a month ( jFalcone et 
|ar ||2010j). While the d ata d id not show a clear peri- 
odicity, [Falcone et al. (2010) determined that periods 



of 35—40 days or > 54 days are possibilities. HESS 
J0632-I-057 has also been associated with the ROSAT 
source IRXS J063258.3-H054857 and the EGRET source 



3EG J0634-I-0521 (Aharonian et al. 2007 Skilton et al 



20091). The source has not yet been detected by t'ermi 
( |Abdo et al.|[20T0l ). A non-detection by VERITAS indi 
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cates a decrease m 7-ray fiu x from the tinie of t he orig- 
inal HESS detection as well ( Acciari et al.|[2009 ). Radio 
observations revealed the source to have a mean flux of 
0.3 mj y at 5 GHz, with va riations on the timescale of a 
month ( Skilton et aH2009 ). The multiwavelength detec- 
tions and the flux variations of this system are consistent 
with the class of 7-ray binaries. 

Currently, six systems are kno wn or proposed 7- ray bi- 
naries: LS 5039, LS I +61 303 ([A bdo et al. |[2C^, PSR 
B125 9-63 ([Aharonian et al.'2005y, Cyg X-'H Albci-t et al ; 
2007| ), Cyg X^ ( .Abdo et al. 2009) , and HESS J 0632+037 
(see above). Optical spectroscopy can provide useful in- 
formation about the stellar properties and orbital pa- 
rameters of these systems, which aid in the analysis of 
observations at other wavelengths. We have undertaken 
a study of the optical properties of HD 259440 to com- 
plement ongoing work in X-ray, 7-ray, and radio to de- 
termine the nature of this system. 

In order to determine the properties of HD 259440, we 
present an extensive collection of red and blue optical 
spectra. Section 2 describes our observations. Section 3 
details our search for radial velocity variations and places 
limits on the orbital period for the system. We address 
the circumstellar disk features in Section 4. In Section 5, 
we discuss our fitting procedure using model spectra to 
determine the stellar properties. Section 6 describes the 
spectral energy distribution (SED) of the Be star. We 
summarize our results in Section 7. 
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2. OBSERVATIONS 

We observed HD 259440 using the Kitt Peak National 
Observatory (KPNO) Coude Feed (CF) telescope with 
the F3KB detector between 2008 October 17 and Novem- 
ber 21. We used both blue and red spectral setups each 
night and generally obtained two spectra of HD 259440 in 
each. The blue spectra were taken using grating B in the 
third order with the 4-96 order sorting filter. They have 
a resolving power R — A/AA ~ 9500, cover a wavelength 
range 4130 — 4570 A, and have exposure times rang- 
ing from 20—30 minutes giving a signal-to-noise, S/N, 
of 50-140. The red spectra were obtained using grating 
B in the second order and the OG550 order sorting fil- 
ter. They have R ~ 12000, span a wavelength range of 
6400 — 7050 A, and have exposure times of 20 minutes 
giving a S/N of 100—250. ThAr comparison lamp spectra 
were taken every 2—3 hours for wavelength calibrations. 
We zero corrected, flatfielded, and wavelength calibrated 
each spectrum using standard procedures in IRAF^ All 
spectra were interpolated to a common wavelengtn grid 
and corrected for heliocentric radial velocity variations. 
The blue spectra were written to a common log wave- 
length scale for radial velocity analysis. The spectra were 
rectified in IRAF using emission line free continuum re- 
gions. 

Two blue optical spectra were obtained at the KPNO 
2.1m telescope using the GoldCam Spectrograph on 2008 
December 12 — 13. Grating G47 was used in the second 
order with a CUSO4 blocking filter and a slit width of 
1.3" to achieve R ^ 2100 — 3100. These observations 
cover a wavelength rang e of 3700-4900 A. The two spec- 
tra have exposure times of 15—20 minutes, resulting in a 
S/N of 125—150. Comparison spectra were taken using 
a HeNeAr lamp to wavelength calibrate the data. The 
spectra were zero corrected, flatfielded, and wavelength 
calibrated using the standard procedures in IRAF. We 
corrected the spectra for heliocentric radial velocity vari- 
ations and rectified them in IRAF using emission line 
free continuum regions. 

We also observed HD 259440 seven times, twice in 2009 
March and five times in 2009 October, with the SOPHIE 
cross-dispersed cchellc spectrograph at the Observatoire 
Haute-Provence (OHP) 1.93m telescope. These observa- 
tions used the 52.65 1 mm~^ R2 echelle grating to cover 
39 orders in the high efficiency mode (i? = 40000) cover- 
ing a wavelength range of 3872 - 6943 A. The exposure 
times were between 45—90 minutes long, yielding a S/N 
of 80—220. Additional spectra were obtained from the 
OHP archive, which were taken between 2007 October 
and 2008 February. These spectra had shorter exposure 
times of 5—15 minutes, and S/N as low as 40. The spec- 
tra were corrected for blaze, flatfielded, wavelength cal- 
ibrated using ThAr lamp spectra taken within 2 hours 
of each observation, and corrected for heliocentric radial 
velocity variations. The observations are summarized in 
Table 111 

3. CONSTRAINTS ON ORBITAL PERIOD 



^ IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



We performed a cross-correlation of the He I A4471 line 
in our CF data with a mean spectrum from that run to 
search for radial velocity, Vu, shifts. Due to the low S/N 
of those spectra, the Vr measurements had a standard 
deviation, (Jj^y, of 20.1 km s^^, and no significant shifts 
in Vft were found. We also did not detect any gradual 
trend of increasing or decreasing Vr. The CF dataset 
rules out orbital periods of 35 d or less. 

The He I A5876 line in the echelle spectra from OHP 
was also examined for Vr variations, of which the results 
are presented in Table [2] These measurements had a 
(Try of 6.7 km s""'^, which is close to the precision that 
we would expect to achieve due to the natural line width 
and line shape. Again, no significant shifts in Vr were 
detected in this set of observations. Due to the signif- 
icant instrumental differences and wavelength coverage 
between the CF, KPNO 2.1m, and OHP runs, we hes- 
itate to combine the V r over thes e comb ined datasets 
to investigate Vr shifts. Crampton ( 1968 ) gives a radial 
velocity measurement of Y.6 km/s; however, the measure- 
ment is flagged for poor quality. Thus we use only the 
sparsely sampled OHP spectra to constrain the orbital 
period further. 

We cannot rule out the possibility that HD 259440 is a 
binary until we determine the probability that the system 
has an orbit with Vr shifts too small for our observations 
to detect. We use a method similar to that outlined in 



Garmany et al. (1980). For a binary system, the mass 
function ot the primary is 



/(to) 



with 



1.0355x10-^ P {1- 



=,2^3/2 



(1) 



Mx 



(2) 



where Mi, is the mass of the optical star in solar units, 
AIx is the mass of the compact object in solar units, Ki, 
is the velocity semi-amplitude of the optical star in km 
s~^, P is the orbital period in days, and e is the orbital 
eccentricity. Assuming an upper limit for if* of 2a rv-, 
we can find the upper limit for the system inclination for 
a given value of P at which we would fail to detect Vr 
variations in our dataset: 

3 fPq{l + q)Hl~e^)y'^'^' 



sini < 9.392 x 10^ 



(3) 

If we assume a random distribution of orbital inclina- 
tions, the probability that Vr is lower than our detection 
limit is 



P{i < iup) 



di 



1 



(4) 



We used the lower aRv obtained from the OHP spec- 
tra and our measured mass of the Be star. Mi, = 16 Mq 
(see Section 5). Secondary masses of 1.5, 2.0, and 2.5 
Mq, consistent with a low mass companion, and eccen- 
tricities of 0.0, 0.2, 0.4, and 0.6 were considered. The 
results are plotted in Figure [T] For a compact compan- 
ion of 1.5 M0, we would be unable to detect Vr shifts 
for an orbital period greater than 50 days. The chance 
of detection improves for higher values of the secondary 
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mass. However, for a 100 day orbit, there is still only 
a 50% probability the the system will have an inclina- 
tion favorable for detecting Vr shifts in our data. The 
chances of detection also decrease with orbital eccentric- 
ity. However, our method docs not take into account that 
to detect the radial velocity shifts for an eccentric orbit, 
we would need to observe the system close to periastron. 
The system would only spend a small fraction of its orbit 
with Vr shifts high enough for us to detect, decreasing 
our chances of observing Vr shifts in systems of higher 
eccentricity. Therefore, we cannot rule out HD 259440 
as a binary system with a long period. 

4. CIRCUMSTELLAR DISK FEATURES 

Our collection of red CF spectra provide an excellent 
tracer of the Be star's disk emission over 35 nights. The 
mean equivalent width of the Ha emission line, H^Ha, is 
—52.3 A, with a standard deviation of 0.7. We find that 
is constant within our estimated 2% error due to 
noise and the continuum placement. We show in Figure 
[2] the Ha line profiles and a gray-scale image of this line 
over the span of the observing run. Since the orbital 
period of HD 259440 is unknown, neither the line profiles 
nor the gray-scale plots are folded by orbital phase, but 
rather they reveal true chronological variations in the line 
profile behavior as a function of heliocentric Julian date 
(HJD). 

We noticed some subtle temporal variations in the 
shape of the emission profile, so we subtracted the mean 
emission line profile to investigate the residuals more 
carefully. The Ha emission residuals, and the corre- 
sponding gray-scales, are shown in Figure[3j These resid- 
ual spectra reveal a partial S-shaped feature over the 35 
nights of observation that suggest a characteristic period 
of ^ 60 days. We attribute these variations to a spiral 
density wave in the circumstellar disk, common among 
Be stars ( [Porter fc Rivinius||2003|) . 

The high V sm i tor this star (see Section 5) implies 
that the disk inclination must be nearly 90°. Therefore 
we view the disk nearly edge-on. Since we do not resolve 
a double peaked profile in Ha, the disk is likely opti- 
cally thick out to several stellar radii (consistent with 
the very high Wna). By contrast, H7 shows a strong 
double-peaked emission profile. The disk also contributes 
a continuum flux to the line of sight which is added to 
the observed stellar spectrum, causing photosphcric lines 
in the rectified spectra to appear too shallow. Bound- 
free emission will contribute a fraction of the stellar flux 
across all wavelengths, while free-free emission will in- 
crease at longer wavelengths. Indeed, we find evidence 
for such flux dilution in our blue spectra, discussed in 
the following section. We also observe an infrared excess 
in the SED for the system, which will be discussed in 
Section 6. 

5. STELLAR PARAMETERS 

Accurately determining the stellar parameters of HD 
259440 is a challenging process due to the prevalence of 
emission features in the spectra. In order to minimize 
the impact of the disk emission on our analysis, we chose 
to focus on the He I and H Balmer lines with A < 4000 
A, where the lines are predominantly in absorption. This 
approach required our analysis to use the lower resolution 
KPNO 2.1m telescope spectra, as these covered a broader 



wavelength range than the KPNO CF data (Figure [1]). 
In order to improve the S/N, we created a mean spec- 
trum from the two observations. We then compared this 
spectrum to Tlust y OSTAR2002 (iLanz fc Huben^ |2003 1 
and BSTAR2006 ( |Lanz fc Hu'beny||2007[ ) model grids, 
corrected for instrumental broadening, to obtain values 
for the effective temperature. Tag, surface gravity, logg, 
and V smi of the optical star. 

The Tlusty BSTAR2006 and OSTAR2002 model 
grids assume a plane-parallel stellar atmosphere. The 
BSTAR2006 models are available with temperatures 
ranging from 15000 K to 30000 K, with a grid spacing 
of 1000 K, and \ogg ranging from 1.75 to 4.75 where 
g is in cm s~^, with a grid spacing of 0.25. The OS- 
TAR2002 models exist for temperatures between 27500 
K and 55000 K with a grid spacing of 2500 K. These 
models are available for log g between 3.0 and 4.75, with 
a grid spacing of 0.25. The OSTAR2002 models as- 
sume a mictroturbulence velocity of 10 km s~^, while the 
BSTAR2006 models assume a microturbulence velocity 
of 2 km s^^. We used models with solar abundances. 

Initially, we tried fitting the mean spectrum of HD 
259440 over the 3750-3950 A range for Tog, log.g, and 
Vsini simultaneously. Tes was varied between 20000 K 
and 40000 K in increments of 1000 K (2500 K) to match 
the spacing of the BSTAR2006 (OSTAR2002) grid. We 
varied logg from 2.50 to 4.75 (3.0 to 4.75) in increments 
of 0.25 according to the BSTAR2006 (OSTAR2002) grid 
spacings. We also increased Vsini from 300 km s~^ to 
500 km s~^ in steps of 50 km s~^. We computed the 
value for each combination of Tcff, log g, and Vsini. 
However, none of the trials produced a satisfactory fit to 
the observed spectrum. The best fits from the two model 
grids produced divergent solutions for Tcff and logy. Fur- 
thermore, a visual inspection comparing the model to the 
observed spectrum revealed the poor quality of these fits 
and they could not be used to constrain a Toff and log g 
for the system. 

The unusually shallow yet broad absorption lines seen 
in the spectra of HD 259440 are indicative of flux dilu- 
tion, in that the continuum emission and scattering from 
the disk has effectively decreased the relative strength of 
absorption lines in the rectified spectra. The total ob- 
served fiux is the sum of the flux from the disk, F^isk, 
and the star, Fstar- Over a small wavelength range, we 
can assume a constant ratio N = Fdisk/ Fstar- Using 
this assumption, we can correct for the disk continuum 
emission by adjusting the rectified spectrum: 



FR,star = Fnil + N) ~ N 



(5) 



where Fr is the rectified observed spectrum and FR ^tar 
is the rectified spectrum of the star corrected for the 
disk continuum emission. Since the relative disk fiux is 
unknown, we adjusted the observed spectrum for values 
of = 0.1 — 0.9 and repeated the fitting procedure at 
each value. 

The absorption lines in the spectrum of HD 259440 are 
fairly broad, so we expect the star to have a high rota- 
tional velocity. At all values of Tgg and logg, a Vsini 
of 500 km s"""^ provided the best agreement with the 
line profiles. Therefore, we chose to fix Vsini at this 
value. We then repeated our fits using the OSTAR2002 
and BSTAR2006 models over the 3750-3950 A region. 
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At each N, we determined the best fit Tcs and log g 
by minimizing the from the model fits. Our quoted 
errors for Tcs and logg represent the la significance lev- 
els from our best fit. For values of N between 0.0 and 
0.4, Tcff and logg are not well constrained, but the fit 
quality and agreement between the BSTAR2006 and OS- 
TAR2002 models improve significantly for higher values 
of N. Therefore, we chose to focus on the results with 
N > 0.5, which are summarized in columns 1—5 of Table 

m 

A contour plot illustrating the Icr, 2cr, and 3(t sig- 
nificance errors for the combined OSTAR2002 and 
BSTAR2006 model fits in the iV = 0.7 case is shown 
in Figure [5j As can be seen from this figure and Table 
(3) the best values of Toff for both the OSTAR2002 and 
BSTAR2006 models lie near the boundaries of the model 
grids. Furthermore, the Icr significance contours are not 
closed at the boundary between the models. Therefore, 
we determined the lower limit of Tog from the la contour 
of the BSTAR2006 models and the upper limit from the 
la contour of the OSTAR2002 models. 

A Tlusty BSTAR2006 model spectrum with Tcff = 
30000 K, log g = 4.0, and Vsini = 500 km s'^ plotted 
against the mean spectrum of HD 259440 with N = 0.7 is 
shown in Figure [61 We generally find good agreement be- 
tween the modeland observed spectrum in the H Balmer 
lines, although the He I A3819 line appears stronger than 
all of our model fits. HD 259440 has been classified as a 
star wit h peculiar abu ndances. It was marked as nitrogen 
rich by Turner ( 1976 1 and may also be a helium-strong 
star. We caniiot distinguish significant differences be- 
tween the model fits with different values of N, hence 
some ambiguity in the physical parameters of this sys- 
tem. 

As a further test of the validity of the fit, we compared 
the model fits to one of the few other photospheric ab- 
sorption lines in our blue spectra, the fiux-corrected He I 
A4471 line. Figure [7] shows good agreement over the core 
of the hne with the BSTAR2006 fit and TV = 0.7 . There 
may be some excess emission in the line wings from the 
higher velocity regions of the disk, narrowing the line 
slightly. Using the OSTAR2002 fits, the flux-corrected 
He I A4471 line is again stronger than the model, a re- 
sult of the difference in microturbulent velocity between 
the two sets of models. Since HD 259440 may be He 
strong, the fit of the He I A4471 line does not help con- 
strain the best fit any further. The quality of the fits 
for the values of N considered in this paper were equally 
good, so we also cannot constrain the fiux contribution 
from the disk any further. 

We determined the mass, M^,, and radius, i?*, of the 
star using our best fit values for Toff and log g and inter- 
polating between the evolutionary tracks of Schallc r^et al.| 
P"992j . The evolutionary tracks were published for non- 
rotating stars with solar metallicity and masses ranging 
from 0.8 to 120 M©. The errors for M^, were calculated 
using the best value of log g and using the la limits for 
Toff, while the errors for i?* were calculated using the best 
fit values for Toff and using the la limits for logg. Our 
measurements of A/^ and i?^ are presented in columns 
6-9 of Table [3| The rapid rotation of HD 259440 may 
result in an underestimation of the mass and introduce 
an additional system atic error into our measurements. 
Ekstrom et al. ( 2008 ) calculated evolutionary tracks for 



rapidly rot ating stars. Compari ng these models to the 
tracks from Schaller et al. ( |1992 ), we conclude that any 
discrepancy in the mass or HD 259440 would be about 5 
Mq. We also calculated the critical velocity, Vcrit, at the 
stellar equator for each N. Since HD 259440 is a rapid 
rotator, the stellar radius at the equator is 1.5-R^ in the 
standard Roche model. Therefore, 



1.5i?/ 



(6) 



The results are listed in column 10 of Table [H These 
results indicate that HD 259440 is rotating close to Vcrit, 
which is plausible for Be stars and thought to contribute 
to the fo rmation of the circumstellar disk (Townsend et 
al] [2004| . 



6. SPECTRAL ENERGY DISTRIBUTION 

We compared the observed SED of HD 259440 with 
the model SEDs from the Tlusty BSTAR2006 and OS- 
TAR2002 models to measure the distance to the star. 
We obtain ed ultraviolet br oadband photometric fiuxes 
|1982| , Johnson UBV magnitudes 
( 1980[ ), Johnson I magnitude from 
7L|MASS JHKs magnitudes from 
!006^, 



from Wesselius et al. 



from Ncckcl & Chini 
^roegc ct al. ( 20071 

'Hkrutskie et al. (2006'), and mid- infrared fluxes from 
|Egan et al., (2003) . The Joh nson magnitudes we re con- 
verted to fluxes according to Bessel et al. (1998), while 
the 2MASS photom etry were converted to fluxes using 



Cohen et al.| ( |2003 ). 

The Tlusty model spectra were binned in 50 A bins to 
eliminate small scale line features and simplify compari- 
son with the observed fl uxes. A reddening of E{B—V) = 
0.85 ( Friedemann||1992 ) and selective extinction of i? = 
3.1 were apphed to the mod el SEDs using the extinction 
model of Fitzpatrick (1999). The observed photometric 
fluxes were corrected tor disk emission using each value 
of N and compared with the best fit OSTAR2002 and 
BSTAR2006 models. 

In our spectral fitting, we assumed that the disk 
flux ratio N is constant over a small wavelength range. 
Therefore, we compared the observed fluxes between 
3000 - 5000 A, corrected for N, to the model SEDs to 
determine the angular diameter, 



obs 



tot. obs 



yR = 



F*, surf y -^^^surf 

(1 + 7V) 



(7) 



where Fi,^obs is the observed flux produced by the star, 
Fi,^surf is the flux at the stellar surface, Ftot,obs is the 
combined flux of the star and disk, and d is the distance 
to the star. Using our derived i?* with the calculated 9fj, 
we find that HD 259440 lies within the distance range 
d = 1.1 — 1.7 kpc. The errors in the distance presented 
in Table[3]give the standard deviation between the model 
SED and the broadband photometric fiux points in the 
3000 - 5000 A range. As the quality of aU our best fits 
were similar for values of N between 0.5 and 0.9, the 
actual error in d is more accurately reflecte d by the range 
of v alues presented. HD 259440 is listed by |Voroshilov et| 
al] (|1985|) as a i nember of the open clu ster NGC 2244. A 



recent study by Bonatto & Bica ( 2009 1 places the cluster 
at 1.6 ± 0.2 kpc from the Sun, a value consistent with 
previous studies and our measurements for HD 259440 
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( ORura fc Ishida||1981[ iPerez et al.||1987| iPark & Sung 

mS|. — — 

We expect to see exorbitant infrared emission in the 
SED of HD 259440 due to the nearly edge-on circumstel- 
lar disk. This infrared excess can clearly be seen in Fig- 
ure [8j which shows the observed broadband photometric 
fluxes with the BSTAR2006 model SED for Toff = 30000 
K, log 5 = 4.0, N = 0.7, and d = 1.30 kpc. At blue op- 
tical wavelengths, the disk contributes ^ 0.7T*. In the 
mid-infrared, the disk continuum emission is ^ 22 times 
the flux of the star. 

7. SUMMARY OF RESULTS 

HD 259440 is a main sequence BOpe star that may be 
helium strong. Due to continuum emission and scattered 
light from the disk contaminating the stellar spectrum, 
our measurements of its physical parameters are some- 
what ambiguous. We find Toff « 30000 K and log g w 4.0. 
From the evolutionary tracks of Schaller et al. (1992), 
our best model fits indicate that the star has a mass of 
about 16 Mq and radius of about 6.6 Rq. We find an 
angular radius of 23.5—25.7 /las and a distance range of 
1.1—1.7 kpc. Due to the non-detection of radial velocity 
shifts over our 35 nights of CF data, we can rule out the 
possible 35—40 d orbital period suggested by Falcone et 
al. (2010). However, longer period orbits would produce 
Vh shifts that would be lost in the noise of our data sets. 
Therefore, we cannot rule out HD 259440 as a binary 
with a period greater than 100 days, which would be 
consistent with the class of non-accreting X-ray binaries. 

We would like to thank Wenjin Huang for allowing 
us to obtain spectra of HD 259440 during his KPNO 
2.1m telescope time in December 2008. We also thank 
Yves Gallant, Bertrand Plezv, and Guillaume Dubus, 
whose observations comprise the OHP archival spectra 
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for support from NASA DPR numbers NNX08AV70G, 
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an institutional grant from Lehigh University. This pub- 
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of Massachusetts and the Infrared Processing and Analy- 
sis Center/California Institute of Technology, funded by 
the National Aeronautics and Space Administration and 
the National Science Foundation. 
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Figure 1. Probability to have missed the binary companion because of a too low inclination angle as a function of the orbital period, for 
eccentricities equal to 0.0, 0.2, 0.4, and 0.6, respectively. In eax;h panel, three different secondary masses have been assumed. The diamonds 
indicate inclinations from 10° to 80°, spaced at intervals of 10°. 



HD 259440: The Proposed Optical Counterpart of HESS J0632+057 




Figure 2. The upper plot shows the Ha hue profile of HD 259440 over our continuous 35 nights of observation with the CF, sorted by 
HJD, and the lower plot shows a gray-scale image of the same line. Since no period is currently Icnown, the "Arbitrary Phase" goes from 
at the start of the observing run to 1.0 at the end of the run. The intensity at each velocity in the gray-scale image is assigned one of 
16 gray levels based on its value between the minimum (bright) and maximum (dark) observed values. The intensity between observed 
spectra is calculated by a linear interpolation between the closest observed phases. 
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Figure 3. The upper panel shows the Ho residual spectra of HD 259440, and the bottom panel shows a grayscale plot of the same 
residuals, in the same format as Fig. [2] The residual spectra reveal variations in the Ha line with a ~ 60 d period that we attribute to a 
spiral density wave in the circumstellar disk. 



HD 259440: The Proposed Optical Counterpart of HESS J0632+057 
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Figure 4. The mean spectra from the KPNO 2.1m (bottom) and CF (top) telescopes are plotted for comparison. The CF spectrum is 
offset by 0.2 for clarity. 
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Figure 5. A contour plot of the Icr, 2cr, and 3<t errors from the spectral fits using N = 0.7 and V sin j = 500 km s~^. The discontinuity 
between the OSTAR2002 and BSTAR2006 models is marked by the vertical dashed line. The square represents the model with the lowest 
value for both model fits. 
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Figure 6. The mean KPNO 2.1m spectrum, adjusted for A'^ = 
model spectrum with T^ff = 30000 K, log g = 4.0, and V s\ni 
shifted by 0.6 for easier comparison. 



0.7, is shown as a solid line. The dashed line shows the Tlusty BSTAR2006 
= 500 km s^^. The residuals from the fit are plotted below the spectrum. 
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Figure 7. The He I A4471 line from the mean KPNO 2.1m spectrum is shown in the same format as Figure [6] 




Figure 8. The Tlusty BSTAR2006 model SED normalized to dji = 24.3 /xas is plotted as a solid line. The diamonds indicate the observed 
broadband photometric fluxes. The triangles represent the estimated stellar flux assuming N = 0.7. 
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Table 1 

Journal of Spectroscopy 



UT Dates 


Telescope 


A Range 


Resolving Power 


Grating/ 


Filter 


Number of 






(A) 


(A/AA) 


Order 




Spectra 


2007 Oct-2008 Feb 


OHP 1.93m 


3872 - 6943 


40000 


R2/50-88 




16 


2008 Oct 17-Nov 21 


KPNO OF 


4130 - 4570 


9500 


B/3 


4-96 


68 


2008 Oct 17-Nov 21 


KPNO OF 


6400 - 7050 


12000 


B/2 


OG550 


62 


2008 Dec 12-13 


KPNO 2.1m 


3700 - 4900 


2100-3100 


G47/2 


CuS04 


2 


2009 Mar 


OHP 1.93m 


3872 - 6913 


10000 


R2/.50 88 




2 


2009 Oct 


Olll' 1.93m 


3872 - ()9 13 


10000 


l-t2/o0 88 
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Table 2 

OHP Vn Measurements 



UT Date Mid-exposure Va 
HJD-2450000 km s'^ 



2007 


Oct 02 


4375, 


,1714 


43, 


,78 


2007 


Oct 22 


4395, 


,1686 


47, 


,55 


2007 


Oct 22 


4395, 


,1744 


46, 


,69 


2007 


Nov 01 


4405, 


,1763 


47, 


,96 


2007 


Nov 02 


4406, 


,1672 


46, 


,48 


2007 


Nov 03 


4407, 


,1905 


44, 


,90 


2007 


Nov 04 


4408, 


,1886 


47, 


,04 


2007 


Nov 05 


4409, 


,1474 


41, 


,74 


2007 


Nov 10 


4414, 


,0348 


47, 


,81 


2007 


Dec 27 


4461, 


,9402 


41, 


,99 


2007 


■Ian 19 


4484, 


,8059 


38, 


,06 


2008 


,lan 21 


4486, 


,8267 


28, 


,88 


2008 


Feb 06 


4502, 


,9462 


30, 


,41 


2008 


Feb 12 


4508, 


,8925 


36, 


,28 


2008 


Feb 19 


4515, 


,8693 


31, 


,02 


2008 


Feb 25 


4521, 


,8527 


27 


,09 


2009 


Mar 10 


4900, 


,8138 


33, 


,11 


2009 


Mar 14 


4904, 


,7721 


45, 


,97 


2009 


Oct 13 


5117, 


,6426 


39, 


,85 


2009 


Oct 15 


5119, 


,6701 


43, 


,57 


2009 


Oct 17 


5121, 


,6604 


41, 


,69 


2009 


Oct 27 


5131, 


,6605 


42, 


,96 


2009 


Oct 28 


5132, 


,6867 


47, 


,04 



Table 3 

Stellar Parameters 



N 






logs 


A log 9 


M* 


AM* 


R* 


AR* 


ycrit 6d 


Ae^ 


d 


Ad 




(K) 


(K) 


(dex) 


(dex) 


(Mq) 


(Mo) 


(Ro) 


(R©) 


(km s~^) (/xas) 


(/ias) 


(kpc) 


(kpc) 



BSTAR2006 Model Fitting Results 

0.5 30,000 tVZ 3.75 tVsl 19.0 9.6 

0.6 29,000 l^iro 3.75 17.6 t^o 9.3 

0.7 30,000 tl'^Z 4.00 +_lfo 16.0 6.6 

0.8 30,000 tlZ 4.00 16.0 tH 6.6 

9 29 000 +^'°°° 4 00 14 5 6 3 

U.» Z»,UUU _6,000 ^-"^ -0.60 '■^■'^ -6.3 -2.9 

OSTAR2002 Model Fitting Results 

0.5 30,000 t^'Zo 4.00 l^o 16-0 6.6 tfg^ 

0.6 30,000 4.00 Igfs 16.0 +JA 6.6 tH 

0.7 30,000 i-^-;;;;;; 4.oo ti^'ri 16.0 tii e.e ti\ 

0.8 27,500 3.75 15.3 8.6 tl% 

0.9 27,500 1^.500 4^00 +0.45 +4.2 +9.4 



501 


25.1 


±1.1 


1.78 


±0.08 


491 


25.2 


±1.2 


1.71 


±0.09 


555 


23.7 


±1.1 


1.30 


±0.06 


555 


23.1 


±1.1 


1.33 


±0.06 


541 


23.4 


±1.2 


1.26 


±0.07 


555 


24.5 


±1.0 


1.25 


±0.05 


555 


23.7 


±0.9 


1.30 


±0.05 


555 


23.0 


±0.9 


l.,34 


±0.05 


476 


24.5 


±1.2 


1.64 


±0.08 


529 


24.1 


±1.3 


1.16 


±0.06 



Note. — The errors in Tgfj and log g were calculated using the OSTAR2002 results to determine the upper limits 
and the BSTAR2006 results to determine the lower limits. 



